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HANDBOOK  OF  SUPERSONIC  AERODYNAMICS 
Volume  4 


Preface 


A  general  preface  to  the  entire  Handbook  of  Supersonic  Aerodynam¬ 
ics  appears  in  Volume  1;  therefore,  the  present  preface  applies  specifical¬ 
ly  to  the  present  issue  of  this  portion  of  Volume  4  only. 

_  ^  A  This  volume,  when  completed,  will  contain  the  following  sections: 
Section  9  -  Mutual  Interference  Phenomena,  Section  10  -  Static  Stability 
Section  11  -  Dynamic  Stability,  and  Section  12  -  Aeroelastic  Phenomena. 
Section  12  only  is  being  issued  at  this  time;  the  remaining  sections  for 
Volume  4  will  be  issued  when  completed. 

Since  the  publication  of  Volumes  1  and  2  the  contents  of  future 
volumes  in  the  Handbook  Series  has  been  changed  in  accordance  with  the  out¬ 
line  set  forth  on  page  iii  of  this  preface  under  the  caption:  "Contents  of 
Future  Volumes  in  the  Handbook  of  Supersonic  Aerodynamics  Series.” 

The  numbering  system  for  Volume  4  is  the  same  as  that  used  in 

Volume  2. 


Agencies  and  individuals  interested  in  the  aeronautical  sciences 
should  feel  free  to  submit  and  recommend  material  for  inclusion  in  the 
Handbook;  full  credit  will  be  given  for  all  such  material  used.  In  the 
selection  of  material  and  the  preparation  of  the  volumes  in  the  Handbook 
Series,  the  Applied  Physics  Laboratory  lays  claim  neither  to  omniscience 
nor  to  infallibility;  therefore,  it  earnestly  solicits  constructive  criti¬ 
cisms  and  suggestions.  Correspondence  relating  to  the  editing  of  the  Hand¬ 
book  Series  should  be  sent  to 


Thomas  F.  Ball,  Supervisor 
Aerodynamics  Handbook  Project 
Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 

Applications  and  communications  concerning  distribution  of  the 
Handbook  Series  should  be  sent  to 


Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 
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SECTION  12 


-  AEROELASTIC  PHENOMENA 


The  following  symbols  are  used  in  the  material  appearing  on  pages 
1200-1  to  1208.2-58  of  Section  12: 


Primary  Symbols 


a 

b 


Ci  >  C, 


"Lh 


'La 


"Mh 


"Ma 


velocity  of  sound  (free  stream),  ft/sec 

semi-chord  length,  ft 

location  of  aileron  hinge  line  measured  from  mid¬ 
chord  point  in  fractions  of  the  semi-chord  (+  aft) 

translational  spring  constant  per  unit  span, 
(lbs/ft)  /  ( ft  span) 

coefficients  of  determinantal  equation 

part  of  supersonic  flutter  aerodynamic  force  co¬ 
efficient  due  to  vertical  displacement  of  the  wing 
quarter-chord  axis  only 

part  of  supersonic  flutter  aerodynamic  force  co¬ 
efficient  due  to  rotational  motion  only 

part  of  supersonic  flutter  aerodynamic  moment  co¬ 
efficient  due  to  vertical  displacement  of  the  wing 
quarter-chord  axis  only 

part  of  supersonic  flutter  aerodynamic  moment  co¬ 
efficient  due  to  rotational  motion  only 


C’ 

^Lh 

CT ,  when 
Lh 

using  the  reduced 

frequency 

of 

the 

aileron 

c* 

'"La 

CLa  when 

using  the  reduced 

frequency 

of 

the 

aileron 

r* 

’"Mh 

CMh  WhGn 

using  the  reduced 

frequency 

of 

the 

aileron 

C’ 

Ma 

CMa  when 

using  the  reduced 

frequency 

of 

the 

aileron 

C” 

Lh 

CLh  when 
forward 

using  the  reduced 
of  the  aileron 

frequency 

of 

the 

wing 

C" 

La 

CT  when 
La 

forward 

using  the  reduced 
of  the  aileron 

frequency 

of 

the 

wing 

C” 

Mh 

CMh  When 
forward 

using  the  reduced 

of  the  aileron 

frequency 

of 

the 

wing 

C" 

cMa 

CMa  when 
forward 

using  the  reduced 
of  the  aileron 

frequency 

of 

the 

wing 

a 


torsional  spring  constant  per  unit  span, 
(  f  t-lbs/rad) /(  f  t  span) 
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Co  torsional  spring  constant  per  unit  span  for  aileron 

p  (  f  t-lbs/rad) /(  f  t  span) 

d  distance  of  elastic  axis  aft  of  quarter-chord  line, 

ft 

E  Young's  modulus  of  elasticity 

Ee  elastic  energy 

E^  kinetic  energy 

F  half  the  rate  of  energy  dissipation 


gh  structural  translational  damping  factor 

gft  structural  torsional  damping  factor 

g structural  torsional  damping  factor  for  aileron 
G  shear  modulus  of  elasticity 

h  displacement  of  wing  quarter-chord  axis  from  the 

neutral  position  (+  downward),  ft;  also  a  general¬ 
ized  displacement 

h*  displacement  of  wing  elastic  axis  from  the  neutral 

position  (+  downward),  ft 

hQ  amplitude  of  h ;  also  generalized  amplitude  of  dis¬ 

placement 

amplitude  of  h* 

i  complex  operator,  ]/ -1 

4 

I  section  moment  of  inertia,  ft 

1^  moment  of  inertia  of  system  about  elastic  axis  per 

unit  span,  lb-f t-sec2/( f t  span) 


1^3  moment  of  inertia  of  aileron  about  hinge  line  per 

unit  span,  lb-f t-sec2/(  f t  span) 

J  effective  section  polar  moment  of  inertia,  ft^ 

k  reduced  frequency,  u)b/V,  non-dimensional 

j^  =  £1  (M2-1)/2M2  j 

k^  reduced  natural  frequency  in  torsion,  ua  b/a 

(,  semi-span,  ft 

L  aerodynamic  force  per  unit  span,  assumed  at  quarter- 

chord  (+  downward,  negative  lift) ^ 


#  — - - 

The  symbol  *L  for  aerodynamic  force,  as  used  in  this  section  of  the 
Handbook,  for  either  primary  or  secondary  concepts,  is  in  the  opposite  direc¬ 
tion  to  that  of  lift  as  customarily  used  in  aerodynamics  and  as  defined  in 
Section  1  of  this  Handbook. 


generalized  aerodynamic  force 

part  of  aerodynamic  force  per  unit  span  (L),  as¬ 
sumed  at  quarter-chord  point,  due  to  various  time 
derivatives  of  vertical  displacement  (h)  of  the 
wing  quarter-chord  axis 

part  of  aerodynamic  force  per  unit  span  (L),  as¬ 
sumed  at  quarter-chord  point,  due  to  rotational 
displacement  of  the  wing 

aerodynamic  force  due  to  aileron  per  unit  span 

mass  of  moving  system  per  unit  span 

mass  of  wing  per  unit  span  ( mj  =  m  in  most  appli¬ 
cations) 

mass  of  aileron  per  unit  span 

Mach  number  (free  stream),  V/a;  also  moment  per  unit 
span  ( +  nose  up) 

generalized  aerodynamic  moment  per  unit  span  about 
elastic  axis 

part  of  aerodynamic  moment  per  unit  span  (M)  about 
the  quarter-chord  axis,  due  to  vertical  displace¬ 
ment  (h)  of  the  wing 

part  of  aerodynamic  moment  per  unit  span  (M)  about 
the  quarter-chord  axis,  due  to  rotational  displace¬ 
ment  of  the  wing 

aerodynamic  moment  about  hinge  line  due  to  the 
aileron 

aerodynamic  moment  per  unit  span,  about  the  elastic 
axis 

/  4 

mechanical  parameter,  I^/Trpb  ,  non-dimensional 

location  of  wing  elastic  axis  measured  from  wing 
mid-chord  point  as  a  fraction  of  the  semi-chord 
(+  aft),  non-dimensional 

mass  unbalance  per  unit  span,  mx^b 

time,  seconds 

air  velocity  (free  stream),  ft/sec 

distance  of  center  of  gravity  chordwise  from  elastic 
axis  as  a  fraction  of  the  semi-chord  (+  aft),  non- 
dimensional 

distance  of  center  of  gravity  of  aileron,  measured 
from  aileron  hinge  line,  in  fraction  of  the  semi¬ 
chord  (+  aft) 

distance  along  span  from  wing  root 

displacement  of  wing  in  rotation  from  the  neutral 
position,  radians/(ft  span),  (+  nose  up) 
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a  displacement  of  wing  in  rotation  from  the  neutral 

J  position,  normalized  in  three-dimensional  case,  per 

unit  span,  radians 

0  angle  of  aileron  with  respect  to  chord  line  of  wing 

(+  trailing  edge  downward) 

A  coefficients  in  the  third  order  stability  equation 

°  (see  Subsection  1207) 

air  density 

Mach  angle  =  arc  sin  l/M  cos^  /j.  =  ( -  1)/M^j 

functions  of  y  defining  the  shapes  of  vibration  modes 

circular  frequency  of  oscillation,  radians/sec 

uncoupled  natural  frequency  in  translation,  y C  /m, 
radians /sec  ' 

u>  uncoupled  natural  frequency  in  torsion, 

radians /sec 

uncoupled  natural  frequency  in  torsion  of  aileron, 
rad  ans /sec 

frequency  parameter,  =  2k/cos^  [i  =  [2M^/(M^  -  1)1  k 


Auxiliary  Symbols 

The  bar  over  a  symbol  (  )  denotes  the  real  component  of  the  com¬ 

plex  quantity  designated  by  the  associated  symbol. 

The  asterisk  (*),  used  as  a  superscript,  denotes  the  imaginary 
component  of  the  complex  quantity  designated  by  the  associated  symbol. 

The  dot  (')  is  used  to  denote  differentiation  with  respect  to 
time,  thus  a  =  da/dt  and  a  =  d^a/dt^. 


P 

P 

*1>  ^2’  ^3 
u> 

“h 
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SECTION  12  -  AEROELASTIC  PHENOMENA 


This  section  of  the  Handbook  of  Supersonic  Aerodynamics  was  pre¬ 
pared  at  the  Applied  Physics  Laboratory  of  The  Johns  Hopkins  University,  with 
the  cooperation  of  the  Bumblebee  Committee  on  Aeroelasticity  and  Structural 
Dynamics.  Members  of  this  committee  were  as  follows: 


M. 

V. 

Barton 

C. 

w. 

Besserer 

H. 

A. 

Cheilek 

M. 

Dubl in 

A. 

H. 

Flax 

H. 

w. 

Pope 

T. 

K. 

Riggs* 

Defense  Research  Laboratory, 

University  of  Texas 

Applied  Physics  Laboratory,  The  Johns 
Hopkins  University  -  Chairman 

Cornell  Aeronautical  Laboratory 

Consolidated  Vultee  Aircraft  Corporation 

Cornell  Aeronautical  Laboratory 

Consolidated  Vultee  Aircraft  Corporation 

Applied  Physics  Laboratory,  The  Johns 
Hopkins  University  -  Secretary 


The  original  draft  of  this  section  was  prepared  for  the  Committee 
by  T.  K.  Riggs  in  accordance  with  the  Committee’s  recommendations  and  sug¬ 
gestions.  The  final  draft  was  prepared  by  C.  N.  Warfield  who  gratefully 
acknowledges  the  helpful  comments  and  suggestions  by  the  members  cf  the  Com¬ 
mittee  and  by  his  colleagues,  F.  K.  Hill,  J.  P.  Kearns,  R.  M.  Mains,  and 
E.  Shotland — and  the  helpful  assistance  of  Mrs.  Corine  Carwile  Bloss  who 
checked  many  of  the  equations  and  the  numerical  results,  computed  the  nu¬ 
merical  example,  and  prepared  the  copy  for  the  final  graphs. 


The  tables  of  flutter  coefficients  which  appear  in  this  section 
were  especially  computed,  under  the  supervision  of  E.  C.  Kennedy,  at  the 
Ordnance  Aerophysics  Laboratory  on  International  Business  Machines  Corpora¬ 
tion  equipment  for  initial  publication  in  this  Handbook. 


1200  Introduction 

1200.1  General  Scope  of  Section 

In  this  section  of  the  Handbook  there  are  presented  certain  tables 
and  graphs  that  may  be  used,  on  the  basis  of  flutter  considerations,  in  the 
design  of  guided  missiles.  In  addition  there  is  included  here  a  brief  treat¬ 
ment  of  certain  theoretical  aspects  of  flutter  in  the  supersonic  regime. 

This  treatment  includes  a  derivation  of  one  of  the  equations  for  flutter  of 
airfoils  in  supersonic  flow,  namely  that  for  torsional  flutter  of  a  two-di¬ 
mensional  (infinite)  wing. 

The  tables  above  referred  to  (Tables  1208.2)  contain  the  real  and 
imaginary  parts  of  the  supersonic  force  and  moment  flutter  coefficients  for 
airfoils.  These  flutter  coefficients  are  equivalent  to  those  originally  de¬ 
fined  by  Borbely  (Reference  12-1). 


* 

Presently  employed  by  Engineering  Research  Associates,  Inc. 
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These  tables  were  computed  by  use  of  a  recursion  formula  that  was 
devised  by  E.  C.  Kennedy,  and  they  are  tabulated  as  a  function  of  a  fre¬ 
quency  parameter  (Q)  for  each  of  several  values  of  Mach  number  (M).  The 
parameter  (Q)  is  related  to  the  reduced  frequency  (k)  and  to  the  Mach  num¬ 
ber  (M)  by  the  equation  Q  =  [2M^/(M^  -  l)]k,  and  a  table  based  on  this  re¬ 
lationship  is  presented  (Table  1208.1).  The  reduced  frequency  is  the  ratio 
between  the  circular  frequency  of  oscillation  ( u>) ,  in  radians  per  second, 
and  the  number  of  times  per  second  that  the  wing,  due  to  its  forward  speed 
(V),  traverses  a  distance  equal  to  its  semi-chord  (b). 

The  tabular  values  for  the  flutter  coefficients  in  the  great  ma¬ 
jority  of  cases  are  believed  to  be  accurate  to  within  one  in  the  last  digit, 
and  in  no  case  is  the  tabulated  value  in  error  by  more  than  two  in  the  last 
digit.  The  Mach  number  range  covered  is  from  1.1  to  12  while  the  value  of 
Q  ranges  from  0.01  to  20.  The  increments  in  both  M  and  Q  are  in  general 
smaller  than  in  existing  similar  tables.  Because  supersonic  flutter  compu¬ 
tations  sometimes  involve  relatively  small  differences  of  coefficients, 
these  coefficients  have  been  computed  and  tabulated  in  most  cases  to  eight 
significant  figures,  although  in  many  applications  three  or  four  digits  will 
suffice. 

Also  included  in  this  section  are  brief  treatments  of  binary  flutter 
(wing  torsion  and  bending  modes)  and  of  ternary  flutter  (wing  torsion,  first- 
and  second-bending  modes,  aileron  and  wing  torsion  and  bending  modes) .  Both 
two-dimensional  (infinite  span)  and  three-dimensional  (finite  span)  airfoils 
are  analyzed.  Brief  discussions  are  given  of  certain  methods  of  solution 
for  the  higher  order  determinantal  equations  that  appear  in  some  of  these 
analyses.  A  brief  mention  of  the  use  of  coupled  and  of  uncoupled  vibration 
modes  in  supersonic  flutter  is  included. 

For  the  purpose  of  familiarizing  the  non-specialist  with  the  tech¬ 
nique  of  flutter  computations,  this  section  includes  a  numerical  example  of 
an  application  of  the  supersonic  flutter  coefficients.  This  example  is  for 
two-dimensional  binary  flutter,  and  is  based  on  the  method  presented  in  the 
Air  Materiel  Center  report  entitled  "Application  of  Three-Dimensional  Flutter 
Theory  to  Aircraft  Structures"  (Reference  12-2). 

In  addition  to  the  list  of  cited  references,  there  is  included  at 
the  end  of  this  subsection  a  bibliography  of  the  more  pertinent  literature 
on  supersonic  flutter. 

The  effects  of  body  motion  and  the  flexibility  of  attachment  of  the 
wing  are  not  discussed  in  this  section  since  these  effects  are  adequately 
covered  in  the  literature  on  svV'Sonic  flutter  (Reference  12-2)  .  Finite  span 
effects,  resulting  in  a  loss  of  lilt  force  at  the  wing  tip,  are  not  taken  in¬ 
to  account;  however,  theoretical  studies  are  available  on  this  subject  (Ref¬ 
erences  12-3,  12-4,  12-5,  12-6  and  12-7).  Empirical  corrections  may  be  used 
to  account  for  tip  effects  with  some  degree  of  reliability. 

The  effect  of  sweepback  on  the  fluctuating  aerodynamic  forces  is 
somewhat  more  complicated  than  the  effect  on  the  static  lift  and  moment  co¬ 
efficients  for  the  same  type  of  wing.  These  sweepback  effects  are  discussed 
in  Reference  12-3.  It  is  possible  to  calculate  the  effects  of  sweepback  on 
the  elastic  properties  of  a  wing  by  the  use  of  approximations,  provided  the 
aspect  ratio  is  sufficiently  high.  Whenever  a  completed  structure  is  availa¬ 
ble  its  elastic  properties  may  be  obtained  from  ground  vibration  tests. 
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1200.2  Basic  Concepts 


An  airframe  at  rest  on  the  ground  in  still  air  will  respond  to  an 
impulse  in  one  of  three  ways.  Depending  upon  the  amount  of  structural  damp¬ 
ing  present  it  will  either  execute  a  series  of  periodic  oscillations  of  di¬ 
minishing  amplitude,  or  return  to  its  initial  state  of  rest  in  the  shortest 
possible  time  (critically  damped),  or  return  more  slowly  to  a  state  of  rest. 

If  the  airframe  at  rest  is  subjected  to  a  sinusoidal  forcing 
function  it  will,  after  passing  through  a  transient  condition,  settle  into 
a  steady-state  vibratory  motion  with  a  frequency  the  same  as  that  of  the 
forcing  function,  and  whose  deflections  and  amplitude  of  vibration  are  de¬ 
termined  by  the  applied  frequency,  as  well  as  by  the  elastic,  inertial,  and 
damping  characteristics  of  the  airframe  structure. 

Since  fluctuating  aerodynamic  forces  result  from  oscillatory  mo¬ 
tions  of  an  airframe,  the  response  of  an  airframe  to  an  impulse  or  sinu¬ 
soidal  forcing  function  will  be  determined  by  these  fluctuating  aerodynamic 
forces  as  well  as  by  the  characteristics  of  the  airframe  structure.  If  the 
phase  relationship  of  the  aerodynamic  forces  is  such  as  to  reinforce  the 
motions  producing  them,  then  a  condition  of  self-sustaining  oscillation  is 
possible.  This  condition  gives  rise  to  what  is  known  as  flutter.  The 
flutter  frequency  is  determined  by  the  flight  Mach  number  as  well  as  by  the 
structural  characteristics  of  the  airframe. 

In  flutter  analyses  computations  are  made  for  the  critical  flutter 
condition  in  which  the  amplitude  of  vibration  tends  to  remain  constant. 

When  the  amplitude  of  vibration  increases  the  condition  is  considered  un¬ 
safe;  when  the  amplitude  decreases  it  is  considered  safe. 

The  boundary  between  the  safe  and  unsafe  flutter  conditions  may  be 
identified  by  investigating  the  equations  of  motion.  An  approximate  meas¬ 
ure  of  the  margin  of  safety  may  be  given  by  the  value  of  the  critical  struc¬ 
tural  damping  factor  computed  for  the  airfoil  structure  with  the  aid  of  the 
herein  tabulated  aerodynamic  flutter  coefficients.  Then  this  value  can  be 
compared  with  the  actual  structural  damping  factor  obtained  experimentally 
by  a  vibration  test,  or  by  estimation  based  on  experience.  Or,  the  degree 
of  safety  from  flutter  may  be  estimated  by  considering  the  distance  between 
the  point  on  a  suitable  chart  describing  the  known  properties  of  the  wing 
and  the  line  on  the  same  chart,  based  on  the  herein  tabulated  flutter  co¬ 
efficients,  which  designates  the  boundary  between  the  "safe"  and  "unsafe" 
regions . 


In  flutter  analyses  the  computations  are  based  on  the  frequency, 
shape  and  phase  relationship  of  certain  vibration  modes  that  are  character¬ 
istic  of  the  structure.  Ideally,  the  principal*  modes  as  they  occur  in 
flight  under  the  aerodynamic  conditions  that  exist  during  critical  flutter 
oscillations  would  be  used  in  these  flutter  analyses.  Theoretically,  in  the 
case  of  three-dimensional  bodies,  there  are  an  infinite  number  of  possible 
vibration  modes.  For  practical  purposes,  however,  the  deformation  of  an  air¬ 
frame  during  a  state  of  critical  flutter  may  bj  assumed  to  be  a  combination 
of  the  deflections  due  to  the  first  two,  three,  or  possibly  four  of  the 
principal  modes  of  vibration — these  principal  modes  correspond  to  the  lower 
frequencies  at  which  the  structure  vibrates  in  resonance.  Approximations 
to  these  desired  modes  may  be  obtained  by  analytical  methods,  or  by  measure¬ 
ments  made  on  the  airframe  while  vibrating  either  at  rest  on  the  ground  in 
still  air,  or  while  in  flight,  or  by  other  experimental  means.  Reference 
12-8  demonstrates  the  feasibility  of  basing  the  analyses  upon  the  actual 
coupled  modes  of  vibration  rather  than  upon  the  fictional  uncoupled  modes. 


Sometimes  referred  to  in  the  literature  on  flutter  as  the  characteristic, 
natural,  or  normal  (coupled)  modes. 
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It  has  been  found  that  flutter  may  occur  in  the  torsional  mode 
without  the  presence  of  a  flexure  component.  This  is  because  at  certain 
frequencies  and  elastic  axis  positions  the  aerodynamic  damping  is  negative, 
that  is,  the  imaginary  component  of  the  aerodynamic  moment  acts  in  phase 
with  the  angular  velocity  so  as  to  accelerate  the  wing  in  rotation  rather 
than  retard  it.  However,  it  has  been  shown  that  such  pure  torsional  flutter 
cannot  occur  at  Mach  numbers  greater  than  1.58  (Reference  12-9)  for  slow 
oscillations,  and  the  limiting  Mach  numbers  for  more  rapid  oscillations  do 
not  differ  much  from  this  slow  oscillation  value. 

If  an  unswept  wing  were  to  oscillate  in  bending  only,  with  no 
rotary  motion,  then  the  aerodynamic  damping  would  always  be  positive,  and 
no  flutter  involving  this  mode  alone  will  occur. 
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1201  Two-Dimensional  Torsional  Flutter 

When  an  airfoil  oscillates  in  a  torsional  mode  only,  various  mo¬ 
ments  about  the  axis  of  rotation  are  involved.  For  a  unit  span  of  the  air¬ 
foil  the  elastic  restoring  moment  will  be  -C^a  (cf.  symbols  list  on  pages 

1200-1  and  1200-3,  and  Figure  1201-1),  and  the  structural  damping  moment  is 
represented  as  a  fraction,  g^ ,  of  the  elastic  restoring  moment,  rotated  in 


Direction  of 


a.  Directions  (The  notation  as  to  directions  is  the  same  as 
that  of  the  NACA  and  the  American  Standards  Association’s 
"Letter  Symbols  for  Aeronautical  Sciences,  Z-10.7,  1950") 


Figure  1201-1  TWO-DIMENSIONAL  WING  NOTATIONS 
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phase  so  as  to  lead  the  latter  by  90  degrees.  The  resultant  of  these  two 
moments  may  be  represented  by  -( 1  +  i  %a)Caa>  where  i  is  the  complex  opera¬ 
tor  Y  -1  .  The  inertial  moment  per  unit  span  is  expressed  by  - 1  '  and  the 

aerodynamic  moment  per  unit  span  about  the  elastic  axis  is  represented  here 
as  M'  .  The  sum  of  these  moments  is  zero,  and  consequently  the  aerodynamic 
moment  may  be  expressed  by 


M'  =  l'aa  +  (1  +  igtf)  Caa  (1201-1) 

Consider  now  the  contribution  to  the  aerodynamic  moment  M'  about 
the  elastic  axis  per  unit  span  due  to  the  rotational  displacement  a  of  the 
wing  from  the  neutral  position.  If  we  let  the  positive  aerodynamic  force 
(that  is,  negative  lift  L  ),  due  to  this  angular  displacement,  act  at  a 
distance  d  forward  of  theaelastic  axis,  and  let  represent  the  aerodynamic 

pitching  moment  about  the  line  passing  through  the  point  of  application  of 
the  aerodynamic  force  Lff,  it  is  obvious  that  such  a  rotational  displacement 
contributes  to  the  moment  about  the  elastic  axis  an  amount  (see  Figure 
1201-2) 


M 


a 


Lad 


L.E. 


T.E 


Figure  1201-2  FORCE  AND  MOMENT  NOTATIONS 


Likewise,  in  view  of  the  effect  of  various  time  derivatives  of  dis¬ 
placement  (h)  of  the  wing  quarter-chord  axis  which  contribute  M^  and  rel¬ 
ative  to  the  quarter-chord,  it  similarly  follows  that  such  a  translatory  dis¬ 
placement  contributes  to  the  moment  about  the  elastic  axis  an  amount 

*h  -  Lh  d 

The  total  aerodynamic  moment  about  the  elastic  axis,  due  to  both  rotational 
and  translatory  motions,  is  therefore 

M’  =  (Ma  ~  La  d)  +  (Mh  ‘  Lh  d) 


( 1201-2) 
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defined  by 


Using  aerodynamic  force  and  moment  flutter  coefficients  that  are 


"Lh 


7rpb3u>3h 


'La 


"Mh 


a 


Ttp 


b3c^o 


M. 


*Pb3<A 


(1201-3) 


'Mo 


a 


4  2 

Tip  b  u>  a 


one  finds  that  Equation  1201-2  becomes 


i  4 

M  =  npour 


C.i  ..A  —  CT  ir  a  +  C,„  tt  -  C 


hd 


"Mcr 


'La  5 


Mh  Td  Lh 


( 1201-4) 


If,  as  is  customary  in  subsonic  flutter  analyses,  we  assume  the  lift  force 
to  act  at  the  quarter-chord  point  then 

d  =  b  (ijr  +  r) 

and  we  find  that  Equation  1201-4  becomes 


=  WJ  [cKa«  -  CLa  +  r)  a.  +  CMh  +  r)  5  -  CLh  *  r)2  * 


(1201-5) 


(Note-  This  equation  for  two-dimensional  flutter  could  have  been 

obtained  directly  from  the  Borbely-Possio  equation  (1203-8) 
by  using  the  relation  h'  =  h  +  Q!d;  cf.  Figure  1201-3.) 


h’ 


Figure  1201-3  DISPLACEMENT  NOTATIONS 
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To  transform  the  motion  parameters  from  the  quarter-chord  axis  to  the  elas¬ 
tic  axis  (see  Figure  1201-3),  let 


h'  =  h  +  ad 


(1201-6) 


For  the  torsional  mode  only  h* 


0;  and  therefore  Equation  1201-6  reduces  to 


(1201-7) 


Equation  1201-5  then  becomes 
» 

M 


a 


CM a  -  CLa  <2  +  r> 


-  C 


Mh 


i  +  r>  +  cLh 


<2  +  r)2 


(1201-8) 


For  harmonic  oscillatory  motion  of  rotation,  we  may  write 


a 


a  e 
o 


itdt 


(1201-9) 


Differentiating  a  (Equation  1201-9)  twice  with  respect  to  time,  and 
substituting  a  and  its  second  time  derivative,  and  Equation  1201-8  into  Equa- 

tion  1201-1,  and  substituting  for  CQ ./ 1^  ,  one  obtains 


i*a)  -!+*?£- 


'Lh 


'Mh 


+  r>2  +  CLa  <*  +  r> 

(2  +  r>]  =  0 


+ 


(1201-10) 


(Note-  This  equation  for  two-dimenSional  torsional  flutter  could 
have  been  obtained  from  the  more  general  determinantal 
equation  for  two-dimensional  binary  flexure-torsion  flutter 
(Equation  1202-9),  by  equating  the  Mg2  +  &22  element  to 

zero,  in  which  Mg2  and  A22  are  defined  by  Equations  1202-7 
and  1202-10,  respectively.) 


For  convenience,  the  real  and  imaginary  parts  of  the  aerodynamic 
coefficient  term  (i.e.,  the  term  included  in  the  brackets)  are  represented 
hereafter  by  Ag2  and  respectively,  whence 


I r  = 

a22 

UMa  " 

CLh 

(J  +  r)2 

+  CT 

La 

(|+  r) 

+  CMh 

<7  *  r) 

and 

A22 

CMa  ' 

CLh 

(2  +  ")2 

+  C* 
La 

<7  * 

+  CJh 

<7  * 

(1201-11) 


The  reason  for  the  use  of  the  subscript  22  will  be  apparent  in  the  subsec¬ 
tion  on  binary  flutter,  1202.  With  this  symbolism,  Equation  1201-10  becomes 


+  lga)  '  1 


(A22  +  iA22) 


0 


( 12U1-12) 
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Equation  1201-12  may  be  written  as  two  equations:  one  including 
only  the  real  terms,  and  the  other  only  the  imaginary  terms.  When  this  is 

done  and  the  substitution  N  =  I^/frpb^  is  made,  the  following  equations  may 
be  obtained: 


( 1201-13) 


g 


a 


N -X 


22 


( 1201-14) 


These  equations  for  two-dimensional  torsional  flutter  may  be  used 
for  a  quick  survey  of  the  flutter  characteristics  of  a  finite  wing  if  one 
first  obtains  an  approximate  spanwise  average  value  for  each  of  the  para¬ 
meters  involved,  e.g.  1^,  b,  r  and  u>a  .  However,  the  use  of  such  spanwise 

average  values  in  the  equations  for  two-dimensional  torsional  flutter  ob¬ 
viously  cannot  be  relied  upon  for  precise  results. 

When  values  of  u)  and  M,  and  therefore  also  of  Ag2  and  for  a 

certain  elastic  axis  location  (r),  are  found  which  satisfy  Equations  1201-13 
and  1201-14,  the  conditions  for  borderline  two-dimensional  torsional  flutter 

are  defined  for  the  conditions  represented  by  the  parameters  I^/irpb^  and 

u^b/a.  The  latter  term,  ^b/a,  is  hereafter  referred  to  as  the  "reduced 

natural  frequency, "  . 

Several  methods  may  be  used  to  obtain  significant  data  from  these 
equations,  two  of  which  are  described  below. 


Method  1.  Computation  of  torsional  damping  factor  g 


a 


(a)  At  each  Mach  number  of  interest,  using  the  mechanical  para¬ 


meter  N  and  the  elastic 
Equation  1201-13,  for  a 
corresponding  values  of 
determined  by 


axis  location  r  of  the  wing,  determine  by  means  of 
series  of  values  of  the  frequency  parameter  O  ,  the 


"a**- 


Then  the  reduced  natural  frequency  k^  can  be 


a 


=  Mk 


(?) 


( 1201-15) 


where  k,  the  reduced  frequency,  is  given  by 

*  ’ 


(b)  Likewise,  by  means  of  Equation  1201-14,  one  can  determine  the 
values  of  ga  corresponding  to  the  same  values  of  fi  that  were  used  in  (a), 

for  the  same  combination  of  values  of  M,  r,  and  N. 

(c)  For  each  value  of  that  was  used  in  parts  (a)  and  ( b)  there 

has  been  obtained  a  pair  of  values  of  g  and  of  k  .  These  pairs  of  values 

or  a 

can  then  be  plotted  as  in  Figures  1201-4a,  b,  c  and  d,  which  represent  four 
combinations  of  fairly  extreme  values  of  r  and  of  N.  Of  course,  figures  of 
this  type  can  be  prepared  for  any  desired  combination  of  values  for  r  and 
N. 
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If  the  borderline  damping  factor  g^  thus  determined  is  negative  or 

less  positive  than  the  actual  structural  torsional  damping  factor  for  the 
structure,  as  determined  by  damped  vibration  test  data,  safety  from  flutter 
is  indicated;  if  it  is  positive  and  greater  than  the  experimental  value,  un¬ 
safe  flutter  is  indicated. 


Method  2.  Computation,  assuming  the  torsional  damping  factor  gQ  is  zero. 

If  one  is  interested  in  determining  only  a  conservative  indication 
of  the  flutter  characteristic  of  the  structure  (that  is,  whether  or  not  the 
structural  parameters  are  such  as  to  indicate  no  flutter  even  if  the  struct¬ 
ural  torsional  damping  factor  g^  is  zero),  then  it  is  necessary  to  determine 

from  Equations  1201-13  and  1201-14  what  combinations  of  the  several  para¬ 
meters  correspond  to  the  conservative  condition  represented  by  g^  =  0.  This 

has  been  computed  for  various  practical  ranges  of  the  several  parameters  and 
the  results  are  given  in  Figures  1201-5.  The  dashed  portions  of  these  curves 
represent  extrapolated  values  only.  In  these  figures  regions  above  the 
curves  are  free  from  flutter,  but  below  these  curves  the  likelihood  of  flutter 
occurring  increases  with  increasing  distances.  For  example,  with  a  struct¬ 
ure  for  which  r=0,  N=20  and  =  0.25,  it  is  evident  that  flutter  is  probable 

only  at  Mach  numbers  between  1.133  and  1.311. 

Other  methods  of  obtaining  and  presenting  results  for  single-de- 
gree-of-f reedom  (torsional)  flutter  are  described  in  References  12-10  and 
12-11. 

The  following  facts  are  important  in  making  a  decision  as  to 
whether  or  not  an  analysis  for  single-degree-of-f reedom  (torsional)  flutter 
is  adequate  in  any  specific  situation: 

(1)  For  elastic  axis  positions  close  to  the  mid-chord,  static 
divergence  (when  second-order  shift  in  aerodynamic  center  location  is  taken 
into  account)  may  be  more  critical  than  torsional  fultter. 

(2)  At  low  supersonic  Mach  numbers  the  flow  may  be  transonic  in 
character,  and  the  applicability  of  linearized  supersonic  aerodynamic  forces 
used  in  these  analyses  would  then  be  in  doubt. 

(3)  For  (u^  / u>a)  <  1,  the  binary  flutter  stability  boundary  will 
usually  be  more  critical  than  these  torsional  ones. 

For  binary  flexure-torsion  flutter  an  approximation  can  be  obtained 
by  the  method  described  in  Subsection  1202;  and  for  actual  finite  wings  more 
reliable  results  can  be  obtained  by  means  of  the  equations  for  three-dimen¬ 
sional  binary  flexure-torsion  flutter  that  are  presented  in  Subsection  1203. 
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Figure  1201-4a  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM  TORSIONAL  FLUTTER; 

g  vs  k  ,  MACH  NUMBER  (M)  INDEPENDENT.  r  =  0  and  N  =  10 

Pi  1  '  ' 


Figure  1201~4b  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM  TORSIONAL  FLUTTER; 

g  vs  k  ,  MACH  NUMBER  (M)  INDEPENDENT.  r  =  0  and  N  =  100 
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Figure  1201-4c  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM  TORSIONAL  FLUTTER; 

g  vs  k  ,  MACH  NUMBER  (M)  INDEPENDENT.  r  =  -1.2  and  N  =  10 
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Figure  1201 -5b  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM 

TORSIONAL  FLUTTER  FOR  ZERO  DAMPING  (g  =0). 
r  =  -0.2  a 
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Figure  1201-5c  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM 
TORSIONAL  FLUTTER  FOR  ZERO  DAMPING  ( g  =  0). 
r  =  -0.4  a 
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Figure  1201-5d  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM 


TORSIONAL  FLUTTER  FOR  ZERO  DAMPING  (g  =  0) . 
r  =  -0.6  “ 
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Figure  1201-5e  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM 
TORSIONAL  FLUTTER  FOR  ZERO  DAMPING  (g  =0). 
r  =  -0.8  a 
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Figure  1201-5f  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM 
TORSIONAL  FLUTTER  FOR  ZERO  DAMPING  (g  =0). 
r  =  -1.0  “ 
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Figure  1201-5g  STABILITY  BOUNDARIES  FOR  SINGLE-DEGREE-OF-FREEDOM 
TORSIONAL  FLUTTER  FOR  ZERO  DAMPING  (g  =0). 
r  =  -1.2  “ 
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1202  Two-Dimensional  Binary  Flexure-Torsion  Flutter 


The  equations  of  motion  for  a  two-dimensional  airfoil  in  flexure 
and  torsion  are  most  easily  derived  (References  12-12  and  12-13)  by  use  of 
the  Lagrangian  equations 


and 


( 1202-1) 


The  quantities  q^  and  q ^  are  the  generalized  coordinates  describ¬ 
ing  the  motion  of  the  system;  they  may  be  considered  as  the  translational 
displacement  h'  of  the  wing  elastic  axis,  and  the  angular  displacement  a, 
respectively,  although  this  choice  is  not  essential.  Thus,  for  harmonic 
oscillatory  motions  we  get: 


h' 


a 


(1202-2) 


The  quantities  L  and  M  are  the  generalized  aerodynamic  force  and  moment 

&  B 

per  unit  span,  respectively. 

The  kinetic  energy  Efe  of  the  system  per  unit  span  can  be  written 

as  the  sum  of  the  translational  and  rotational  energies  about  an  axis  through 
the  center  of  gravity,  as  follows, 

2 

Ek  1  j  "  [h'  +  *ab  «]  +  I  [xa  "  m  (xt*b)2]  «2  (1202-3) 

Expanding,  substituting  S  for  the  mass  unbalance  quantity  mx^b,  and  also 
writing  the  equations  for  the  elastic  energy  Ee  and  half  the  rate  of  energy 
dissipation  F  per  unit  span,  one  obtains: 

Ek  =  |  (mh’2  +  2Sh’a  +  i'a®2) 

Ee  =  \  (Chh'2  +  C ac?)  (1202-4) 

F  =  1  (5L^  h’2  +  g*C«  a2) 

By  introducing  the  generalized  coordinates  q^  and  qg  (Equations  1202-2)  into 

these  energy  equations,  taking  derivatives,  and  then  substituting  into  the 
Lagrangian  equations  of  motion  (Equations  1202-1),  we  have: 
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2  .»  icut 
u)  mhoe 


2„  iu)t 
SaQe 


chh;e^‘  *  ig„chh;e 


iurt: 


-  L_  =  0 


2. .  iu>t  2„.  <  iujt  _  iu)t 

w  Shoe  +  Caa oe 


ig  C  a  e 
&a  a  o 


iiut 


-  M„  =  0 


( 1202-5) 


The  generalized  force  and  moment  per  unit  span  on  a  two-dimensional 
wing  about  the  elastic  axis  (see  Equations  1203-10)  are: 


_  T  t  _  .3  2  iu>t 

-  L  =  -  TTpb  to  e 

(An 

11 

o 

b 

+  A12 

.  w  ,4  2  iu>  t 

-  M  -  -  Trpb  to  e 

(A21 

h’ 

o 

b~ 

+  A22 

“o) 


(1202-6) 


where  (see  Equations  1203-9) 


A11  “  CLh 


A 


CtK  U  +  r)  -  C 


12  Lh  v2 


La 


(1202-7) 


A21  CLh  ^2  +  CMh 


'22 


CMa  ~  CLh  v2 


4  *  ■'I2  *  <CLa  *  CMh)  (1  *  r) 


By  combining  Equations  1202-5  and  1202-6,  and  rearranging  (since 


=  KC h/m,  u>a  =  Vca  /  1^  ,  and  S  =  mx^b),  we  have: 


In  order  for  a  solution  to  exist,  the  determinant  of  Equations 


1202-8  must  vanish.  That  is, 

M11  + 

A11 

M12 

+  A12 

M21  + 

A21 

M22 

+  A22 

0 


(1202-9) 
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where 


'll 


igh> 


1 


M12  =  M21 


mx 


7rpb 


M22 


+  iga)  - 


1 


( 1202-10) 


Some  methods  of  solving  the  determinantal  equation  for  two-dimen¬ 
sional  binary  flutter  will  be  covered  in  Subsection  1204.  The  determinantal 
equations  of  motion  derived  here  (Equations  1202-7,  1202-9  and  1202-10)  are 
identical  to  those  presented  in  Reference  12-14. 
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1203-1 


1203  Three-Dimensional  Binary  Flexure-Torsion  Flutter 

Let  the  quantities,  h'  and  a,  describing  the  motion  oi  the  three- 
dimensional  (finite)  wing  referred  to  the  elastic  axis  be  defined  by  (cf. 
Equations  1202-2) 


h'  =  qj  =  V;eiu;t 

“  =  *2  q2  =  VoelWt 


(1203-1) 


where  and  4>  ^  are  functions  of  the  spanwise  position,  y.  The  quantities 

q^  and  qg  are  generalized  coordinates;  they  may  be  considered  respectively 

2S  .the  displacement  of,  and  rotation  at,  the  tip  of  the  wing,  although  in 
any  specific  case  some  other  quantity  may  be  more  convenient 

The  kinetic  energy  in  such  a  system  may  be  found  from  the  span- 
wise  integration  (cf.  Equation  1202-4) 


r  i  0 

1  l 

T  fmh,2dy  + 

0 

2^  Sh’ady  +  |  I 'abP dy J 

The  elastic  energy  Ee  in  such  a  system  is 


E  =  A 

e  Z 


| EI  (01)2<,i’ +  JGJ  (ff)2,Jy] 


(1203-Zb) 


One-half  the  rate  of  energy  dissipation  is 
F  = 


■*[-£/ «(#)** -S'!"  $)**] 


(1203-2c) 


Since  h’  and  a  have  been  defined  in  Equations  1203-1,  derivatives 
necessary  for  substitution  in  Equations  1203-2  may  be  formed.  After  sub¬ 
stitution  we  have: 


Ek  = 


E  =  w 
e  2 


r  i  i  i  f 

W  xldy  +  2  |  S^^qjqgdy  +  J  ^  dyJ 

s[  fEI“i2  (^)  “*  +  | 


( 1203-3) 


*  i [  -  51  i2 (^r)  «y-gzr  |  ^2  (P-)  “s'] 


The  Lagrangian  equations  of  motion  for  such  a  system  of  two  degrees 
of  freedom  are  (cf.  Equations  1202-1): 


d 


/dEk\ 

(  r~)  + 

'  9Qi  / 


dE 


9QJ 


e  +  9F 


9<li 


-Lg  =  ° 


(1203-4) 
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d 

3T 


+ 


^  '“g 

aq2 


0 


where  and  are  the  generalized  aerodynamic  force  and  moment  per  unit 

span  acting  on  the  wing,  referred  to  the  generalized  coordinates  q^  and  q2, 

respectively.  The  former  will  be  more  fully  defined  in  Equations  1203-7 
and  1203-8,  respectively. 


Taking  the  necessary  partial  derivatives  of  the  energy  equations 
(1203-3)  and  substituting  into  the  Lagrangian  equations  (1203-4),  we  have: 


2  icot  .  i 
-  u)  e  h 


t 

o  |  m<*> i 


2  ,  2  iwt 

dy  —  id  e 


l  1  d2  </>  2 

«0  J  s  «i  hdy  +  hoeiwt  | EI  (-j-y-)  dy 


1  rd2*i  2 


*  ^yut  |  si  (^)  <>*-*■,= 


=  0 


(1203-5) 


2  iwt  .  . 
—  U)  0  n 

o 


J  S  *2dy  -  a)2 eiajt  ao  f  I^«22dy+  ffQeiaJt  f  GJ  (hT^)  dy 
0  8  0 


+  igaaoeiu,t  I  GJ  (^)2<*  -  *g  “  0 


(1203-6) 


Borbely's  and  Possio's  equations  for  the  lift  and  moment  on  a  unit 
span  of  two-dimensional  wing  oscillating  in  flexure  and  torsion  are  derived 
in  References  12-1  and  12-15,  respectively,  and  are  reproduced  in  Reference 
12-14.  Using  the  coefficients  defined  by  Equations  1201-3,  the  force  and 
moment  about  the  elastic  axis  may  be  written,  respectively: 


L  = 


-  i3  2  iwti  ~o 

-  itpb  u)  e  {_cLhF-  + 


M'  =  -  7rpb4w2elwt 


<Z  +  r>  CLh  -  CL«  ] 

['  CHh  +  +  r)  CLh]  B2  +  [  ‘  c 

-  CLh  <Z  *  r)2  *  CLa  (Z  +  r)  +  °Mh  <Z  +  r)]  ”o 


Mo- 


(1203-7) 


(1203-8) 


(Note-  This  equation  for  M'  is  derived  independently  in  Subsection 
1201;  see  Equation  1201-5  and  the  note  that  follows  it.) 
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For  convenience,  let 
A, 


11 

CLh 

12 

=  CLh 

<1 

+  r)  - 

CLa 

21 

=  CLh 

+  r)  - 

CMh 

22 

CMa 

- 

CLh 

+  r) 

(1203-9) 


2  +  C  f1 
+  CLa  2 


.  /I 
"Mh  '7 


ien,  for  two-dimensional  wings, 


T  *  .3  2  iwt 

L  =  -  up  b  u)  e 


*  .4  2  iut 

M  =  -  npb  u)  e 


(AH  F°  +  a12«0) 


-  hQ  . 

,A21  +  A22“o) 


( 1203-10) 


For  three-dimensional  wings,  taking  into  account  the  spanwise  variations  of 
displacement  (cf.  Equations  1202-2  and  1203-1),  we  have 


L  -  -  up 


b3a>2e 


iojt 


(An 


^  +  A 

+  Ai2 


M’  =  -  7rpb4w2eiu,t  (a21  +  A22  *2  aQ) 


h  ao) 


(1203-11) 


By  the  principle  of  virtual  work,  and  by  use  of  Equations  1203-1 
and  1203-11,  the  generalized  moments  and  forces  may  then  be  expressed  as 
follows : 


T  2  iwt 

L  =  -  npw  e 


Mg  =  -  Jrpoj 


ho  J  b2An^dy  +  aQ  J  b3Aj2  ^dy] 

0  0 

?el"t  [ho  |  *1  +  “o  J  b^A22  *2  <*] 


(1203-12) 


0  0 
These  may  be  substituted  into  Equations  1203-5  and  1203-6,  respec¬ 
tively,  to  obtain  the  equations  of  motion,  thus: 


(MH  +  A^)  hQ  +  (M12  +  A 1 2 ^  ao  ~  0 

(l*2i  +  A21^  +  +  A 99)^  =  ® 


( 1203-13) 


22  22'  o 


A  necessary  condition  for  the  existence  of  a  solution  of  these 
equations  is 


t 

11 

+  Aj  j 

M12  +  Ai2 

21 

+  Ail 

M22  +  A22 

=  0 


( 1203-14) 
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where 


M* 

m22 


A22 


up  |  b3A214»1<>2dy 


b4A 


22*2 


dy 


(1203-15) 


In  general,  for  three-dimensional  wings,  each  factor  in  every  one 
of  the  foregoing  integrands  is  a  function  of  its  spanwise  location,  for  vari¬ 
ous  reasons  as  indicated  below: 


Wing  Characteristic  Determining  the 
_ Spanwise  Function _ 


Quantities  So 
Determined 


Mass  distribution 
Material 

Cross-section  form 
Planform 

Planform  and  elastic  axis  location 
Mode  shape  in  flexure 
Mode  shape  in  torsion 


■’  S>  la 
E,  G 

I,  J 

b 

All»  A12’  A21’  A22 
*1 


Further,  it  is  seen  that  the  quantities  M^,  M^,  M21’  and  M22  are 

functions  of  the  mechanical  parameters  and  frequency,  but  not  of  the  flight 
conditions.  However,  the  aerodynamic  terms  Ajj,  A^,  Ag^  and  Agg,  are  func¬ 
tions  of  Mach  number  and  the  location  of  the  elastic  axis  relative  to  the 
mid-chord  line,  as  well  as  of  the  frequency  and  certain  mechanical  parameters. 


For  special  cases,  the  above  equations  may  be  simplified  to  a 
large  extent;  for  instance,  a  uniform  rectangular  cantilever  wing  would  ena¬ 
ble  the  computer  to  remove  all  terms  other  than  0-  and  <j>n  from  the  inte¬ 
grands.  1 


Several  methods  of  solving  the  determinantal  equations  (e.g.  Equa¬ 
tions  1202-9  and  1203-14)  are  possible.  A  method  based  on  that  of  the  U.  S. 
Air  Force  Air  Materiel  Command  (Reference  12-2)  is  presented  as  an  example  in 
Subsection  1204. 
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1204  Applications  of  Determinantal  Equation  for  Two-Dimensional  Binary 

Flutter 


1204.0  Discussion 

The  determinantal  equation  for  two-dimensional  binary  flutter  (cf. 
Equation  1202-9)  is 


M11  +  A11 


M21  +  A21 


M12  +  A12 


M22  +  A22 


=  0 


w‘?ere,  (cf.  Equations  1202-7  and  1202-10): 

[(30* 
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“21 
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igh>  -  1 


mx 


ffpb 
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22 


3W 


(1  +  ig^)  -  1 


A11 = 

CLh 

A12  = 

CLh  (2 

+  r>  *  CLa 

A21  = 

CLh 

+  r>  ~  CMh 

A22  = 

CMa  ' 

‘  CLh  +  r) 

'airly 

simple 

solutions  to 

( 1204.0-1) 


(1204.0-2) 


'La 


'Mh 


equation  have  been  obtained,  and  one  of  these  is  outlined  in  the  following 
subsection . 


1204.1  Materiel  Center  Method  (Peferences  12-2  and  12-16) 


Let 


and 


Then 


( 1204.1-1) 


( 1204.1-2) 
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The  determinantal  equation  may  therefore  be  written 


+  C{  A  +  C2  =  0 

where  Cj  and  C2  are  complex  constants. 


( 1204.1-3] 


The  two  complex  roots  of  this  quadratic  equation  are  given  by 


-Ci  -  4C2 


(1204.1-4) 


By  complex  algebra  it  is  readily  shown  that 


/C1  ~  4C2 


+  V  2  (cos  |  +  i 


(1204.1-5) 


where  £  =  (Cj  -  4Cg) 

r)  *  ( C12  -  4C2)* 

Q  =  arc  a*'  ~ 

Hence  we  may  write  the  real  and  complex  parts  of  the  two  roots  of  Equation 
1204.1-4  as  follows: 


2A1  =  “  C1  + 


2A2  “  C1 


f2  +  v2  cos  | 


'  9 9  Q 

^  +  r\l  cos  ^ 


(1204.1-6) 


2Al  =  ‘  C1  +  V  +  sin  i 


2  9 


2A2  »  -  Cj  -  y  l;2  +  T]2  sin  | 

By  the  definitions  of  Equations  1204.1-1  it  is  apparent  that 

(¥)2  >  A, 

(1204.1-7) 


&2)' 


and  since  co  must  be  assumed  in  order  for  values  of  the  aerodynamic  coeffi¬ 
cients  to  be  chosen,  then  the  values  of  u>  .  are  determined  for  this  value  of 

’  a 

co  and  for  the  simultaneously  assumed  value  of  Mach  number  M. 
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It  is  also  apparent  by  the  definitions  of  Equations 
the  damping  coefficients  are: 


* 


1204.1-1  that 


(1204.1-8) 


Thus,  the  procedure  for  determining  the  stability  of  the  wing  at  a 
given  Mach  number  consists  of: 


(a)  Assuming  a  series  of  values  for  the  reduced  frequency  k,  there¬ 
by  determining  the  values  of  the  frequency  parameter  and  of  the  aerody¬ 
namic  coefficients  which  will  be  used  in  the  determinantal  equation;  then 
using  these  coefficients  in  solving  the  determinantal  equation  for  the  natural 
frequency  in  torsion  and  for  the  damping  factor. 

( b)  Plotting  these  computed  damping  factors  against  some  con¬ 
venient  parameter  such  as  or  to^b/a. 


(c)  Determining  experimentally,  or  estimating  from  experience,  the 
actual  damping  factors  of  the  wing;  and  plotting  this  factor  on  the  graph 
referred  to  in  (b). 

If,  at  a  particular  value  of  Mach  number  and  natural  frequency  cu  , 
the  actual  damping  factor  of  the  wing  is  greater  than  the  computed  value  a 
(i.e.,  if  the  point  representing  the  experimental  value  lies  above  the  curve 
representing  the  computed  values)  then  freedom  from  flutter  is  indicated. 


1204.11  Numerical  Example  by  the  Materiel  Center  Method 


Let  the  following  values  be  assumed  to  define  the  characteristics 
of  a  two-dimensional  wing  that  is  to  be  examined  for  binary  flutter: 

— -  »■  =  100.0 

7rpb^ 

I' 

- 2L  =  16.67 

7rpb 


CtV 

7T  =  0.700  (1204.11-1) 

a 

r  =0 


These  values,  when  substituted  in  the  M-terms  (Equations  1204.0-2) 
of  the  determinantal  equation  give: 

Mn  =  100.0  (0.4900  A  -1) 

m12  =  M21  =  0 


( 1204.11-2) 
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M22  =  16.67  ( A-l) 

Let  the  flight  Mach  number  (M)  of  interest  be  1.4;  and  let  the  fre¬ 
quency  range  of  interest  be  defined  by  a  range  from  0.2  to  0.7  for  the  fre¬ 
quency  parameter  d  .  For  this  immediate  part  of  the  numerical  example  the 
value  0.4  is  chosen  for  the  latter  quantity. 

That  is 


M  =1.4 
fl  =  0.4 


(1204.11-3) 


These  two  values  determine  the  aerodynamic  coefficients  (as  tabu¬ 
lated  in  Table  1208.2)  to  be: 

CLh  =-1.31345  -  i  12.999891 


CL(J  =-132.93679  -  i  6.776163 
CMh  = -1-08389  -  i  6.367874 


(1204.11-4) 


CMa  =-fi5- 34705  +  i  3.340791 


For  r  =  0  and  for  these  coefficients,  the  A-terms  (Equation 
1204.0-2)  of  the  determinantal  equation  become: 


A11 

A12 

A21 

A22 


=  1.313  +  i 
=  132.280  - 
=  0.427  -  i 
=-1.335  +  i 


13.000 
i  13.276 
0.132 
0.113 


(1204.11-5) 


Substituting 
the  A-terms  (Equations 
get: 


these  values  for  the  M-terms  (Equation  1204.11-2)  and 
1204.11-5)  into  the  determinantal  equation  1204.0-1,  we 


49.00  A-  98.69  +  i  13.00 
0.4272  -  i  0.1321 


132.28  -  i  13.28 

16.67  A  -  18.00  +  i  0.1133 


This  equation  when  expanded  and  simplified  gives 
A2  +  (-3.094  +  i  0.2721)  A+  (2.106  -  i  0.2719)  =  0 


=  0 

(1204.11-6) 

(1204.11-7) 


By  comparison  of  this  equation  with  Equation  1204.1-3  it  is  appar¬ 
ent  that  the  complex  constants  are: 


Cj  =  -3.094  +  i  0.2721 


C 


2 


2.106  -  i  0.2719 


(1204.11-8) 
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The  quantities  that  appear  in  the  roots  of  the  determinantal  equa¬ 
tion  can  be  calculated  by  Equations  1204.1-5  as  follows: 


Cl2 

= 

9.4994  -  i  1.6838 

4 

C2 

= 

8.4258  -  i  1.0876 

Ci2  -  4 

C2 

= 

1.0735  -  i  0.5962 

i 

= 

1.0735 

(1204.11-9) 

As 

V 

= 

-0.5962 

■  s 

= 

1.1081 

e 

= 

-29°  2.74’ 

The  real  and 
equation  are  therefore, 

imaginary  parts  of  the  two  roots 
by  use  of  Equations  1204.1-6: 

of  the  quadratic 

Al 

= 

2.083 

A2 

= 

1.011 

(1204.11-10) 

* 

A1 

= 

-0.2750 

* 

a2 

= 

0.0029 

By  Equations  1204.1-7  it  is  apparent  that  the  natural  frequencies 
of  the  wings  {ioa)  in  relation  to  the  circular  frequency  (u>)  corresponding  to 

the  specified  Mach  number  and  frequency  parameter  are  given  by: 


=  1.443 

U) 

(1204.11-11) 

war2 

=  1.005 

0) 

A  convenient  non-dimensional  parameter  for  the  natural  frequency  is 
u)a  b/a;  this  can  be  derived  from  the  foregoing  ratio  by  the  identity 

wab'  .wa  »b  V 
a  -  ~  w  ‘  V  *  a 


where  wb/v  ( the  reduced  frequency  k)  is  related  to  Mach  number  M  and  fre¬ 
quency  parameter  Q  as  indicated  in  the  list  of  symbols  and  in  Table  1208.1. 


and 


For  this  numerical  example  we  therefore  find  that 


ub 

IT 

V 

a 


0.09796 

1.4 


(1204.11-12) 


NAVORD  REPORT  1488  (Volume  4)  HANDBOOK  OF  SUPERSONIC  AERODYNAMICS 

1204.11-4  Aeroelastic  Phenomena  1  January  1952 


Therefore, 


“a  lb 


wa2b 


0.1980 


=  0.1379 


( 1204.11-13) 


By  Equations  1204.1-8  it  is  apparent  that  the  damping  coefficients 
corresponding  to  the  two  roots  of  the  flutter  equation  are: 


gj  =  -0.1320 
g2  =  +0.0029 


(1204.11-14) 


Similar  computations  of  o^b/a  (the  reduced  natural  frequency  kff), 

and  of  g,  have  been  computed  for  A  =  0.2,  0.25,  0.3,  0.5,  0.6,  0.7,  and  1.0 
and  then  all  of  these  values  have  been  plotted  (g  vs  )  in  Figure  1204.11-1, 

for  Mach  number  1.4.  In  an  actual  investigation  of  the  flutter  characteris¬ 
tics  of  a  wing  similar  computations  and  graphs  would  be  computed  for  each  of 
several  other  Mach  numbers. 

The  actual  value  of  the  quantity  k^  for  the  sample  wing  may  be  de¬ 
termined  by  experiment  or  estimated  from  experience.  In  the  former  case  the 
natural  frequency  in  torsion  of  the  wing  (wa  in  radians  per  second)  would 

be  measured,  and  also  an  average  or  effective  semi-chord  length  of  the  wing 
would  be  determined.  In  addition,  for  each  altitude  of  interest  a  value  for 
the  velocity  of  sound  would  be  determined  corresponding  to  the  ambient  tem¬ 
perature  and  composition  of  the  air  at  that  level. 


Likewise  the  actual  value  of  the  damping  factor  of  the  wing  in 
torsion  would  be  determined  by  measuring  the  rate  of  decay  of  a  damped  tor¬ 
sional  vibration  of  the  wing  structure,  or  by  measuring  the  power  required 
to  sustain  such  a  vibration  at  constant  amplitude  -  or  an  estimate  could  be 
made  of  the  torsional  damping  factor  from  past  experiences.  A  similar  de¬ 
termination  would  be  made  of  the  flexural  damping  factor  of  the  wing  struct¬ 
ure,  and  both  of  these  damping  factors  would  be  used  in  selecting  a  suitable 
common  damping  factor  for  the  wing  being  considered. 

The  point  representing  the  value  of  the  damping  factor  (g)  and  of 
the  non-dimensional  parameter  for  the  reduced  natural  frequency  ( )  of  the 

wing  at  a  given  altitude  would  then  be  plotted  on  the  previously  computed 
graphs  such  as  represented  in  Figure  1204.11-1  for  each  Mach  number  of  in¬ 
terest.  If  the  point  for  the  experimental  quantities  lies  above  both  curves 
representing  the  two  roots  of  the  equation  it  is  concluded  that  flutter  is 
improbable.  For  example,  if  the  quantity  kff  for  the  wing  at  sea  level  is 

0.2527  and  the  smaller  of  the  two  damping  factors  is  0.0032  it  is  seen  that 
the  point  representing  this  wing  on  the  graph  of  Figure  1204.11-1,  for  M  = 
1.4  lies  above  both  curves  and  therefore  the  wing  appears  to  be  free  from 
flutter  at  this  Mach  number.  Similar  spotting  of  the  experimental  values  on 
the  graphs  for  other  Mach  numbers  would  be  made  to  determine  the  possibility 
of  flutter  occurring  at  each  of  these  Mach  numbers. 


Figure  1204.11-1  ROOTS  OF  EQUATIONS  DETERMINING  STABILITY  BOUNDARY  FOR  BINARY 

FLEXURE-TORSION  FLUTTER.  MATERIEL  CENTER  METHOD.  M  =  1.4 
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1205  Three-Dimensional  Ternary  Flexure-Flexure-Torsion  Flutter 

(References  12-12  and  12-13) 

In  many  cases  of  three-dimensional  systems  it  will  be  found  that 
the  natural  frequency  in  bending  in  the  second  mode  may  be  nearly  equal  to 
the  natural  frequency  in  torsion.  If  this  is  found  to  be  true,  then  it  may 
be  expected  that  the  second  bending  mode  will  affect  the  flutter  character¬ 
istics.  In  order  to  include  the  effects  of  the  additional  bending  mode, 
let : 


h  =  hl  +  h2 
hj  =  <^(y)q1(t) 

h2  =  ^(y)  q2(  *) 

a  =  03(y)q3(t) 


*2h2o 

Voe 


i<ot 

e 

iu  t 
e 

iu>  t 


(1205-1) 


The  process  of  determining  the  kinetic  and  elastic  energies  of  the 
system,  taking  appropriate  derivatives  and  substituting  in  the  Lagrangian 
equations  of  motion,  can  be  followed  as  in  Subsection  1203.  If  this  is  done, 
the  condition  that  the  equations  of  motion  have  a  solution  will  be 


where 


M11 

+ 

a"  m"  + 

A11  M12 

a" 

a12 

"is  + 

A" 

A13 

*21 

+ 

a”  m"  + 

21  m22 

a" 

a22 

-is  + 

a” 

a23 

=  0 

(1205-2) 

-31 

+ 

A"  If”  4. 

31  *32 

a" 

a32 

-33  + 

a" 

a33 

-i'l 

= 

U2[(^)2<‘ 

+  ighl> 

-  ljdy 

1 

Mj2  =  Mg j  =  -  j*  m^^gdy  =  0  (by  orthogonality) 

1 

m13  =  M31  =  "  J  s^i^3dy 

o  (1205-3) 

*22  =  |  m*22  [(~w)  < 1  +  igh2)  '  J]  dy 

i 

*23  =  M32  =  '  |  S^2^3dy 

*33  =  |IflV[&)2(1  +  V-1]dy 
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Aii  =  nP  J  b2Ai i^i2 dy 


l12  =  A21  =  b  All*l*2  dy 


* 

a13  =  71  p  f  ^12*1*2  dy 

°  l 

A22  =  np  J  h2A j i  $2  dy 


1205-4) 


A23  ~  77(3  \  t>3A12<^2  ^3  dy 

0 

l 

A31  np  |  b  A21 ^1^3  dy 
o 

l 

A32  =  Trpj  b3A21^^3  dy 


A33  ~  71 P  [  b^A22^3^dy 


The  values  of  the  unprimed  Ajj,  A^2,  A21»  A22  are  tbe  same  as  in 

Subsection  1202  (Equation  1202-7).  The  method  of  solving  Equation  1205-2 
will  be  discussed  in  Subsection  1207.  An  application  of  this  general  method 
to  subsonic  flutter  is  given  in  References  12-12  and  12-13. 
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1206-1 


The  determinantal  equation  for  two-dimensional  ternary  bending- 
torsion-aileron  flutter  may  be  written,  corresponding  to  that  for  binary 
flutter  (Equation  1202-9),  as 


11 

+ 

A11 

M12 

+ 

A12 

M13 

+  A13 

21 

+ 

A21 

M22 

+ 

A22 

M23 

CO 

CM 

< 

+ 

=  0 

(1206-1) 

31 

+ 

A31 

M32 

+ 

A32 

M33 

+  A33 

where  Mjj . ^22  and  ^11 . ^22  are  exactly  as  defined 

in  Subsection  1202.  These  are  repeated  here  for  convenience.  In  addition, 
the  forces  and  moments  about  the  elastic  axis  due  to  the  motion  of  the 
aileron,  and  the  moments  about  the  aileron  hinge  line  also  are  given  here. 
Thus, 


[£): 


(1  +  igh)  -  1 


“12  =MZ1="^ 


“13  1  *31  '  '  ^ 

M23  =  M32  =  "  (c  ~  r> 

7Tp  b  TTpb 


(1206-2) 


=  . 

^  [(?)  +  V  -  >] 

"33 

The  aerodynamic  coefficients  not  involving  the  aileron 
given  in  Equations  1202-7,  that  is: 

are  identi- 

A11  ~ 

CLh 

A12  = 

CLh  <Z  *  r>  -  CLa 

( 1206-3) 

> 

CO 

II 

CLh  <5  +  r>  -  C«h 

A22  = 

CMa  -  CLh  *  r^2  +  (CLa  +  CMh^ 

(j  +  r) 
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The  aerodynamic  terms  involving  the  aileron  are: 

’  '  <TT>  (I  CLh  *  CLar) 

■  '  (i2£)  [cila  *  CLh  (H“l  +  t)  +  CLa  (2  f^-T  +  s)  +  2  c»b] 

*  CLh  <2  *  <=>  -  -  (L^)4  (I  CLh  -  CHh) 

1  -  c„„  -  CLh  <2  +  r>  <2  +  c)  ♦  CLa  (J  +  c)  -  C,h  (>  *  r) 


(1206-4) 


[-  - 1  (2  - 1) +  C  (*  H4  -2)  +I CL>] 

a33  =  '  ( 1  2""C *  (cila  +  ?  CLb  +  2  CLa  +  2  c»h) 

All  of  the  aerodynamic  flutter  coefficients  (i.e.,  all  of  the  C,  Cr 
and  C"  coefficients)  are  obtained  from  Table  1208-2,  in  which  values  of  the 
coefficients  are  tabulated  with  Mach  number  (M)  and  the  frequency  parameter 
(  fl  )  as  independent  parameters,  where  the  latter  is  a  function  of  M,  V,  uj, 
and  b  (see  the  symbols  list).  In  the  case  of  the  C-coef f icients,  b  is  the 
semi-chord  of  the  entire  wing;  for  the  C* -coefficients,  b  is  the  semi-chord 
of  the  aileron;  and  for  the  C"-coeff icients,  b  is  the  semi-chord  of  that 
portion  of  the  wing  forward  of  the  aileron.  For  any  given  wing-aileron  com¬ 
bination  it  is  assumed  for  flutter  analyses  that  the  circular  frequency  u) 
is  the  same  for  all  primed  or  umprimed  C-coef f icients . 

It  should  be  noted  that  if  the  aileron  flutter  alone  (with  no  wing- 
torsion  or  bending)  is  being  investigated,  the  two  families  of  curves  in 
Figures  1201-4  and  1201-5  apply,  if  the  aileron  is  assumed  to  be  hinged  at 
the  leading  edge  (i.e.,  r  =  -1.0). 
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1207  Solution  of  Higher  Order  (above  second  order)  Determinantal 

Flutter  Equations 


If,  in  the  ternary  flutter  determinantal  equations  of  motion  (e.g., 
Equations  1205-2  and  1206-1),  it  is  assumed  that  the  frequencies  bear  a 
fixed  ratio  to  each  other,  and  that  structural  damping  factors  are  equal,  we 
may  write: 


Z 


2 


g 


Sh  =  8 


Of 


0 


(1207-1) 


A  =  Z  ( 1  +  ig) 


It  is  then  found  that  the  ternary  determinantal  equations  may  be 
put  in  the  form  of  a  third  degree  polynomial  such  as 


AqA3  +  AXA2  +  A2A+  A3  -  0  (1207-2) 

Since,  in  supersonic  flutter  analyses  it  is  necessary  to  solve  the 
determinantal  equation  for  each  Mach  number  of  interest,  it  is  obvious  that 
considerable  computational  work  is  required.  Three  methods  of  solving  these 
higher-order  flutter  equations  (including  quadric  as  well  as  cubic  equations) 
have  been  investigated  by  Ruggiero  and  recorded  in  Reference  12-17. 

As  an  alternative  to  solving  the  cubic  equation,  one  may  assume 
that  the  bending  and  aileron  frequencies  are  fixed  quantities  instead  of  being 
in  fixed  ratios  with  the  torsional  frequency.  Then,  on  expanding  the  deter¬ 
minant,  the  stability  equation  will  be  linear  in  A  and  the  torsional  fre¬ 
quency  may  be  found  directly.  After  plotting  u>a  and  gfl  versus  l/k  or  some 

other  parameter,  it  will  be  found  that  at  some  value  of  k  the  calculated  cu 

will  be  the  same  as  the  actual  natural  frequency.  Thus,  the  torsional  damp¬ 
ing  factor  found  at  that  value  of  k  will  determine  the  stability  of  the 
system. 


Other  modifications  of  the  method  may  be  made,  for  instance:  (1) 
assume  the  aileron  natural  frequency  known,  and  the  value  of  kft  known,  and 

then  solve  the  resulting  quadratic  in  A  ;  (2)  assume  that  the  damping  is 
zero,  the  aileron  natural  frequency  known,  and  then  solve  for  Z  and  k^  . 

These  methods  may  also  be  applied  in  principle  to  the  binary  equations  dis¬ 
cussed  in  Subsection  1204. 
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Reduced  Frequency  ( k) ;  Mach  Number  ( M)  and  Frequency 
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